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Abstract. In order to produce appreciable amount of primordial black holes (PBHs), the 
square amplitude of curvature perturbation must take a large value of O(0.01), namely, seven 
digits larger than the value observed by cosmic microwave background radiation (CMB) on 
large scales. Such a large fluctuation can be achieved by violating the slow-roll (SR) condition 
within a short duration. The best known of such possibilities is the ultraslow-roll (USR) 
inflation. We calculate the power spectrum of curvature perturbation in a simple single-field 
inflation model which evolves through the SR-USR-SR regimes so that both large-amplitude 
small-scale fluctuation for PBH formation and small-amplitude large-scale fluctuation as 
observed by CMB are realized. We further calculate the bispectrum and one-loop correction 
to the power spectrum induced by the third-order action of curvature perturbation as the 
beginning of precision cosmology on small scales. As a result, we show that single-field 
inflation model realizing PBH formation can be constrained by the quantum correction. 


Contents 
1 Introduction 1 


2 Two-point functions 2 
2.1 Slow-roll approximation 4 
2.2 Ultraslow-roll inflation 8 
2.3 Constant-roll inflation 9 


2.4 Ultraslow-roll to slow-roll transition 10 
2.4.1 Smooth transition 10 

2.4.2 Sharp transition 11 

3 Three-point functions 12 
3.1 Slow-roll approximation 13 
3.2 Ultraslow-roll inflation 14 
3.3 Ultraslow-roll to slow-roll transition 15 
3.3.1 Sharp transition 15 

3.3.2 Smooth transition 16 

4 Generation of large fluctuation 17 
4.1 Power spectrum 17 
4.2 Bispectrum 20 
4.3 Beyond tree-level contribution 21 

5 Discussion and summary 23 
5.1 Discussion 23 
5.2 Summary 25 


1 Introduction 


PBHs are formed when an overdense region with a magnitude of density fluctuation close 
to unity enters the Hubble horizon during radiation domination [1-4]. The most popular 
formation mechanism of such a large perturbation on small scales relevant to PBHs is quan- 
tum effects during cosmic inflation in the very early universe [5-7] (see [8] for a review). In 
this chapter, we focus on single-field inflation models that can produce large fluctuations on 
small scales relevant to PBH formation with particular attention to our recent finding [9, 10] 
that higher-order quantum effects can restrict the feasibility of such models, even though 
they are observationally possible because there are no stringent observational constraints on 
primordial fluctuations on small scales [11-19] unlike on large scales where observations of 
CMB tightly constrain the fluctuation spectrum [20-22]. 

The simplest inflation model consistent with CMB observation is canonical SR inflation, 
where the scalar field which drives the inflation, called inflaton, is described by a canonical 
kinetic term and a potential. At a certain range, the potential has to be slightly tilted to 
realize the SR condition, so it can explain CMB observations. Beyond such a range, one can 
construct a theory that violates the SR condition, which usually results in the amplification 
of fluctuation on small scales [23, 24] due to rapid slowdown of the inflaton’s velocity. The 


simplest example of a violation of the SR condition is the USR condition [25-29], which can 
be generalized to the constant-roll condition [30-32]. It occurs when the inflaton enters an 
extremely flat region in the potential [33], where the first derivative of the potential vanishes. 
Peaks of the enhanced fluctuation can collapse into PBH after such a fluctuation reenters the 
horizon. An example of potential that realizes the formation of PBHs while satisfying the 
observational constraint on large scales is shown in Fig. 1. 

In this chapter, we explain how violation of the SR approximation can lead to ampli- 
fication of primordial fluctuation on certain scales. Although we choose to focus on USR 
inflation, the method presented in this chapter can be easily generalized to other mechanism 
that yields rapid slowdown of the inflaton’s velocity. In Sec. 2, we introduce two-point func- 
tions or the power spectrum of the fluctuations in SR inflation, USR inflation, and USR to 
SR transition. In Sec. 3, we introduce three-point functions of fluctuations in the same setup 
as in Sec. 2. In Sec. 4, we work out an inflation model with SR-USR-SR transitions. In this 
model, fluctuations are amplified on scales corresponding to the time of the USR phase. In 
Sec. 5, we discuss other PBH formation mechanisms within a single-field inflation framework 
and summarize our review. 


2 Two-point functions 


The action of canonical inflation is given by 


S=5 | av= [MAR - (0,0)? -2V(6)], (2.1) 


where Mp is reduced Planck scale, g = det guv, gu, and R are metric tensor and its Ricci 
scalar. Consider a spatially flat, homogeneous, and isotropic background, 


ds? = —dt? + a?(t)dx? = a?(r)(—dr? + dx’), (2) 


where 7 is conformal time. Equations of motion for the scale factor a(t) and the homogeneous 
part of the inflaton ¢(t) are the Friedmann equations 


H? = : 


des 
= = H =-—> 2: 
saz (58+ V0), i= -zf (2.3) 
p p 
with H = a/a being the Hubble parameter, and the Klein-Gordon equation 

.. av 
3H —=0. 2.4 
Ga BO a (2.4) 


Here, a dot denotes time derivative. 
During inflation, the evolution of the Hubble parameter is parameterized by the slow-roll 
(SR) parameters. The SR parameters are defined as 
én H 


En+1 = en H’ q = H2’ (2.5) 


where n is a positive integer. In case the Hubble parameter H is almost a constant, we have 
€1 < 1. This condition is called the quasi-de Sitter (quasi-dS) approximation. Thus, the 
scale factor can be approximated as 


1 
ax -— xe, (2.6) 
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Figure 1. Schematic picture of the inflaton potential realizing PBH formation from USR condition. 
When the inflaton is around ¢cmp, scales probed by CMB observations leave the horizon and it is in 
the SR regime. It enters an extremely flat region at t = t, undergoing an USR period. It enters the 
SR period again at t = te until ena, the end of inflation. 


Another stronger condition is called the SR approximation, which is defined as |e,| « 1 for 
every n. Therefore, the SR approximation implies a quasi-dS approximation, while the con- 
verse statement is not true. Violation of the SR approximation does not imply a breakdown 
of the quasi-dS approximation. 

Small perturbation from the homogeneous part, (t), of the inflaton (x,t) and metric 
can be expressed as 


p(z, t) = p(t) + 64(z, t), 
ds? = -N?°dt? + qiz (dz* + N‘dt)(da? + NIdt), (2.7) 


where 7; is the three-dimensional metric on slices of constant t, N is the lapse function, and 
N’ is the shift vector. We choose comoving gauge condition 


6b(x,t) =0, yj(@,t) = aP (He) S, (2.8) 


where ¢(a#,t) is comoving curvature perturbation. N and N’ are obtained by solving the 
constraint equations. In this review, we do not consider the tensor perturbation. 
Expanding the action (2.1) up to the second order of the curvature perturbation yields 


Ss [¢] = Ma, f dt dz a®e, G = Oc] . (2.9) 


Introducing the Mukhanov-Sasaki variable by v = zMpi¢, where z = ay2e1, the action 
becomes canonically normalized 


II 
Slol = 3 far ate |e’? - Oar +e], (2.10) 


where a prime denotes derivative with respect to r. We can see that the quantity 2”/z 
behaves as an effective time-dependent mass squared. In momentum space, quantization is 


performed by promoting the Mukhanov-Sasaki variable to an operator 
dn (7) = Mpiz(r)Ce(7) = vn (7) Ge + ož (r)ât ps (2.11) 
where mode function v(T) satisfies the Mukhanov-Sasaki equation 
ut + G = =) Up = 0. (2.12) 


In terms of the SR parameters, the effective mass term can be written as 


1 


3 1 1 1 
| = (aH? (2 ey + 58 5 E12 + 12 jee) : (2.13) 


zZ 


From (2.10) we find that 7” = v’ is momentum conjugate to v, so canonical quantization 
imposes a commutator 


[on(T), te (7)] = iħô(k — k”), (2.14) 
and the operators satisfy [Gz, al] = (27)°6(k — k’) under the normalization condition 
Ue UE — vv% = th. (2.15) 


The two-point functions of curvature perturbation and power spectrum can be written as 
(Gie(7) Caer (T)) = (277)?5°(& + k’) (Ce (7) C0 (7))) , (2.16) 
k3 
A?(k, 7) = a MGa(ICw(M)), (2.17) 


the bracket (---) = (0|---|0) denotes the vacuum expectation value (VEV), and A?(k,7) is 
the power spectrum multiplied by the phase space density. 
2.1 Slow-roll approximation 


With SR approximation, the effective mass becomes 


n 


Z = (aH) |2 — & + Žo + O(e*)| . (2.18) 


Up to first-order correction in SR approximation, aH as function of 7 reads 


1 


H = —-—_——_.. 2.19 
a r(1 = €1) ( ) 
Then, the effective mass can be expressed in terms of the conformal time as 
ð 1 3 Lf. l 
> a2 (2 +361 + se) = ( = i) ’ (2.20) 
where v is defined as 3 3 i 
=L Sya = 65: 2.21 
n= Tim ee (2.21) 


At early time, T + —oo, the Mukhanov-Sasaki equation approaches 


vt + ku, = 0, (2.22) 


which is the equation of a simple harmonic oscillator. Requiring minimum energy state and 
normalization (2.15), solution of the Mukhanov-Sasaki equation at early time is given by 


ja 
= IR 


This minimum energy state corresponds to the initial vacuum state called Bunch-Davies 
vacuum, defined as Gx |0} = 0. 
At general time, solution of the Mukhanov-Sasaki equation is given by 


ee, (2.23) 


uR(T > 


ue(t) = VZT [er HS (kr) + 2HP (—kr)], (2.24) 


where cı and c2 are coefficients to be fixed to a boundary condition. Given the asymptotic 
limit of Hankel functions 


im HO) (a) = Jew [+i (z = = = *)| : (2.25) 


we can match the general solution to the early time limit (2.23) that leads to 


[T 4x 
G = T and c2 = 0. (2.26) 


Then, the mode function of Mukhanov-Sasaki variable becomes 


u(t) = af FeO) THP kr) (2.27) 


To obtain the mode function of the curvature perturbation, we need to know the explicit 
function of z(7). Starting from 


z' 1 
—=aH|1+- 2.28 
= =aHt (1+ 52), (2.28) 


we can substitute (2.19) to obtain 


JŽ J= 1 ) 
= 1 +é + =e 
Z T 2 


2(r) T o 


Tx 


(2.29) 


where 7, is an arbitrary time reference and z, = z(7,). Choosing 7, as horizon crossing time 
k = a(T4)H (T4) so 


(2.30) 


va) » (2.31) 


then the explicit function of z(T) is 


z(T) = (Fv) [-kr (1 — ay ss) (2.32) 


Then, the mode function of curvature perturbation reads 


o R(T) H af eitseo iT rola 
S(T) = TR ain) = (ome), Vi pr) tatia) iiD (FH (kr). (2.33) 


With another asymptotic limit of Hankel function 


lim HP (e) = -tro (Z) ”, (2.34) 


x—>0 T 


the curvature perturbation at the end of inflation T = 7(—> 0), approaches 


B iH 1,3 
Ce (70) = ( Mn an) or °T(3/2). (2.35) 


Also, the power spectrum reads 


A3(k, To) = (sta) (2.36) 


21F2 
87 Mie 


with a weak k-dependence due to the horizon crossing condition manifested in the spectral 

tilt dlog A2(k, 70) 
og A; (E, To 

k l= 

ms(k, 70) dlog k 


Another way to derive power spectrum in SR approximation is by looking at exact de 
Sitter case. In exact de Sitter, the effective mass is 


= ( 2€1 €2)* : (2.37) 


2. (2.38) 


vk(T) = a (1 = x) a, (2.39) 


When we proceed to derive the curvature perturbation, we run into a trouble because in exact 
de Sitter we find e1 = 0. However, as a first-order approximation of the SR parameters, we 
can simply consider the finite value of €;. Also, the Hubble and SR parameters have a weak 
time dependence. 

Numerical calculation of ¢;,(7) shows that evaluating the proportionality coefficient 1/z 
to vk(T) at the horizon crossing time T = 7, yields an excellent fit. So the curvature pertur- 
bation (2.33) can be approximated as 


— iH Le —ikT 
(T) = (z wa) 0 + ikT)e f (2.40) 


Substituting 7 > 0 leads to the power spectrum (2.36). 
The power spectrum (2.36) can be expressed as a function of potential. We can rewrite 
the Friedmann equation as 
V = HMA — 61). (2.41) 


Taking derivative with respect to ¢@ yields 
dV 1 


V = — = —H°’ M (6 =2 ; 2.42 
s450 A pl(6 — 261 + €2) (2.42) 
At leading order in the SR approximation, the first SR parameter can be written as 
M2 V 2 
pl K 
a a r 2.43 
TEO (2.43) 


Therefore, €; has an implicit time dependence because ¢ depends on time. 

From (2.36), we can see that the power spectrum at the wavenumber p corresponds 
to é evaluated at r = —1/p. This implies that observation of the power spectrum at 
the wavenumber p constrains V.4(¢(r)) at rT = —1/p. To solve the horizon problem, the 
fluctuations in the CMB must be around 60 e-folds before the end of inflation. Moreover, 
fluctuations are observed within the finite range of wavenumber 0.005Mpc~! < p < 0.2Mpc7! 
by CMB. This means that CMB observations constrain only a finite region of the potential 
around 60 e-folds before the end of inflation. Beyond such a region, there are much looser 
observational constraints on the potential. Therefore, it is possible to have an inflation model 
in which the SR approximation is violated after the horizon-crossing time of the CMB scales 
and before the end of inflation. 

Comparison with observations can only be made if and only if we are convinced that 
quantum fluctuations generated during inflation behave as classical statistical fluctuations to 
leave observable traces in our universe. The following intuitive argument can be developed 
for the classicalization of curvature fluctuations. In superhorizon regime, k < aH, the mode 


function behaves as 
iH k? 
= 1+ 0| = ]], 2.44 
&lT) (a wa), | + (m) (2.44) 


which means ¢ (7) = —C;(7) to the leading order. So we find ĈklT) = Cy (T) (âk — a!) and its 
conjugate momentum ae (T) = M2.27(7)C,(T) (Ge — ât ,) have the same operator dependence, 
and apparently commute with each other, behaving classically [34]. This is why long-wave 
quantum fluctuations behave as if they are classical statistical fluctuations to provide the 
origin of large-scale structures and CMB anisotropy. The more precise statement, however, 
is that the quantum commutator 


[d.(r), 2g (r)| = ihò(k — k’) (2.45) 
always holds, and it is the commutator 
X z B ih 
[éo(r), Ca(r)| = aay? 9 (2.46) 


that decreases exponentially during the standard slow-roll inflation. 

Observation of the spectral tilt n, = 0.9649 + 0.0042 on the CMB scale implies that 
the SR parameters €1,|€2| at 7 = —1/p have a value O(0.01), based on (2.37). If the first 
SR parameter decreases rapidly at a later time tT >> —1/p, the power spectrum is amplified 
at wavenumber k >> p. This can be achieved by violating the SR approximation at a later 
time. In the next subsections, we will show some inflation models that violate the SR 
approximation. 


2.2 Ultraslow-roll inflation 


In this subsection, we consider a regime where the potential is constant, called ultra-low- 
roll inflation (USR) [25-29]. Because dV/d¢ = 0, the Klein-Gordon equation (2.4) becomes 


$ = —3H¢, so ¢ x a~. This makes «€; strongly time-dependent and extremely small as 
A -6 
€1 = 2M2H? xa 5 (2.47) 


and the second SR parameter becomes 


€1 $ 
=Á = 24 +2- x —6, 2.48 
€2 ZH € + ZH (2.48) 


which breaks SR approximation. The effective mass can be written as 


z" 2 3 La 2 

A = (aH) 2+ 562+ 7+ O(«1) = = (2.49) 
which is the same as in the exact de Sitter case. Thus, solution of the Mukhanov-Sasaki 
equation with Bunch-Davies initial vacuum is the same as (2.39). Then, the curvature per- 
turbation is given by 


G(T) 1+ ikre". (2.50) 


so 
7 2MpivV/ ke1(7) 


Here, €; has a strong dependence on time, which can be expressed as 


e1(7) = e1 (Tx) ( i, i: (2.51) 


where 7, is the horizon-crossing time of a characteristic wavenumber k,. Therefore, the 
curvature perturbation becomes 


(T) - ( ae ) (1 + ikr)e~**7 (2.52) 
T)= : : 
i 2M) ke, (Tx) a(Tx) 
At the end of inflation, the power spectrum is 
H? a(t) \° 
A?(k, To) = l 2.53 
A , To) 8r? M3 (Tx) (>) ( ) 


To solve the horizon problem, CMB-scale fluctuations must cross the horizon at Nomp © 
60 e-folds before the end of inflation. It means that power spectrum of CMB fluctuations 
can be written as , 1? om 
= CMB 
As (p, To) = 8r? M2 e1 (Te) e . (2.54) 
To explain the observed A?(p, 7) = 2.1 x 107°, for a reasonable value of €1(7,), an extremely 
small value of H/M,) is required. Because H/M,) characterize the energy scale of inflation, an 
extremely small value of H/M,) can reach the energy scale of the Big Bang nucleosynthesis. 
Hence, inflation with only USR period cannot have the required e-folds number to solve the 
horizon problem. 


2.3 Constant-roll inflation 


We can generalize the USR condition to €2 = constant [30-32]. By taking reference at time 
Tą, the first SR parameter can be written as 


e1(7) = €1(%) ( ate) i" ; (2.55) 


a(T) 
Rewriting the second SR parameter as €2 = —2(3+ a), the effective mass can be expressed 
as 
2 1 1f 1 
where v is defined as , 
v=ja+ ol. (2.57) 


Solution of the Mukhanov-Sasaki equation with Bunch-Davies initial condition is the same 
as (2.27). The power spectrum is given by 

H? 1 
87? M2 e1 (T0) 27 


3—2v 
A? (k, T) = 27 T (v)]? (=) (2.58) 


We can read that the spectral tilt is ns — 1 = 3 — 2v. Therefore, there are two possible 


solutions for a, given by 

3 1 
Q+ = =) ae z = Ns). (2.59) 
However, we have to choose a solution that does not grow outside the horizon. Recall 


the second-order action of ¢ that reads 


sO] = M3, I dr Be ares [(C)? — (a¢)?). (2.60) 
The equation of motion is given by 


(aĉa) 


acı 


Ch + Ge + k?Ce = 0 (2.61) 
At superhorizon scale, we can neglect the last term in the equation of motion. Hence, 
evolution of the curvature perturbation outside the horizon is 


culo) > Ar + Bef 2 4 


; een) (2.62) 


where Ay and By, are time-independent coefficients. The second term determines whether 
Ck(To) grows or decays. Substituting (2.55), the integral becomes 


TO dr 69+3 
i, Pa mo a (2.63) 


It shows that curvature perturbation grows outside the horizon if e2 < —3 or equivalently 
a > —3/2. It was pointed out in [24] that the curvature perturbation grows outside the 
horizon if 3 — €; + €2 < 0 is satisfied. Clearly, USR condition satisfies this condition since 
gQ= —6. 


= sharp 
= smooth 
= sharp (PBH) 


Figure 2. Sketch of €2(7) for sharp and smooth transition from USR to SR period are shown by 
the gray and blue curve, respectively. Brown curve shows sharp transition from SR-USR-SR period, 
which will be discussed in Sec. 4. 


From observation of spectral tilt ns = 0.9649 + 0.0042 on CMB scale, we obtain a4 ~ 0 
and a_ % —3. Solution a, = 0 corresponds to e2 ~ —6, which approximately equal to 
USR inflation. As explained in the previous subsection, although it yields an almost scale- 
invariant power spectrum, this solution cannot explain the required number of e-folds of 
inflation because the curvature perturbation grows outside the horizon. Therefore, we have 
to choose a_ as the solution because it corresponds to decaying curvature perturbation. With 
this solution, constant-roll inflation can explain the observed CMB fluctuations. 


2.4 Ultraslow-roll to slow-roll transition 


In this subsection, we consider an inflation model in which the second parameter SR evolves 
from the USR to the SR period at conformal time Ts. For T < Ts, the inflaton is in the USR 
period. At time T > Ts, the inflaton undergoes a transition from the USR period to the SR 
period. In the next part, we consider two different evolutions of the second SR parameter: 
smooth and sharp transitions. The former means that the second SR parameter evolves as a 
continuous function of time, while in the latter it evolves as a step function of time. We follow 
the parameterization of €2(7) in [35], which is also used in many literatures. The sketches of 
€2(T) for both cases are shown by the gray and blue curves in Fig. 2. 


2.4.1 Smooth transition 


At time T > Ts, consider evolution of the first and second SR parameter as 


Te\ S 2 

3/(=) (s—3)+s+3 7 \st3 

iai | (=) : (2.64) 
2s(s + 3) 


(ZY 6-3) +s+3 


(2.65) 


where s = y9 -— 12ny ~ 3 -— 2ny and ny —> 0. The parameter ny corresponds to the 
second derivative of the potential at the transition point. In this model, the second SR 
parameter evolves from USR period €2(T;) = —6 to SR period €2(7) = —2ny. With this 


— 10- 


parametrization, the effective mass becomes 


z! 1 1 1 
gR (2 — 3nv) = = (v = z) , (2.66) 


where 
v? = — — 38ny. (2.67) 


Because ny — 0, solution of the Mukhanov-Sasaki equation with Bunch-Davies initial vac- 
uum is the same as (2.39). Then, the curvature perturbation is given by 


G(T) 1+ ikre T, (2.68) 


7 iH ( 
o 2 Mpv k3e1(7T) 


and power spectrum at the end of inflation becomes 


A?(k, To) = SE (2.69) 
, 82? M2 e (T0) 
2.4.2 Sharp transition 
Consider evolution of the second SR, parameter as a step function of time 
e2(T) = —Ace2 [1 — O(T — Ts)|, (2.70) 


where Aen = 6. At time T < Ts, the inflaton is in the USR. period with the Bunch-Davies 
initial vacuum. In this period, the curvature perturbation and its time derivative are given 


by 
iH 


—tkr . 
T) = —————e 1+ ikr 
Ce(7) Ma J Ba) ( ) 
1H Ts 3 —ik Š 
= e "(1+ 7tkr), 2.71 
2Mpi y k’e1 (Ts) ( T ) ( ) ( ) 
I iH Ts\ 3 ai ae ) 

= e k*r — — — 3ik |. 2.72 
L o (2.72) 


At time T > Ts, the inflaton is in the SR period. In this period, the curvature perturbation 
and its time derivative are given by 
iH —ikt . ikT . 
Chr) = [Are (1+ ikr) — Bge“ (1 -— ikr)| : (2.73) 
k3e1 (Ts) 


Ce(T) ker (Are - Bye'* ) (2.74) 


E 2Mpiy/ k3€1 (75) 
where A, and By are general coefficients to be determined from boundary conditions [36-43]. 
This time, we have to consider both positive and negative frequency solutions because the 
boundary condition at 7 = 7, does not correspond to Bunch-Davies vacuum. Solutions of 
the coefficients by requiring continuity of (T) and ¢;,(7) at transition T = Ts are 


3(1 + kr?) 


Ap=12: ee (2.75) 
3(1 + ikTs)? vikr, 


— 11- 


At the end of inflation, the power spectrum becomes 


H? 


Apa til. 2.77 
MZe) k k a 


k3 
2 _ 2 
As (K, 70) = zal) = 
For a very long-wavelength modes, |kTs| < 1, the power spectrum approaches 


ie Ate os 
[rs] Tal” Sr Mge (Ts) 


3 Three-point functions 


So far, we have considered the second-order action of the curvature perturbation that leads 
to two-point functions or the power spectrum at the tree level in the language of quantum 
field theory. Within the cosmological perturbation theory framework, higher-order correla- 
tion functions resulting from higher-order expansion of the action are expected to be small. 
However, because we are interested in PBH formation, where the curvature perturbation is 
large on certain small scales, higher-order correlation functions can be quite large. 

In this section, we discuss the third-order action that yields three-point functions of the 
curvature perturbation. Expanding (2.1) to third-order of ¢ yields the interaction action [44] 


SOI] = Soule] + Sele] sm fat ate 200) (F) (3.1) 
The bulk interaction Spux[¢] reads 
Sbuik[¢] = m f at daa ‘lee c+ m O:C)?C — 2e1Å8;¢3ix 
-zaë F a¢(0:ð;? T zaace] , (3.2) 


where y = €0~2C. Sp[¢] is the surface action given by the total derivatives [45, 46]. The 
last term is interaction proportional to the equation of motion in the lowest order, 


($), =< (€,a3¢) — €a07¢. (3.3) 


The function f(¢) is explicitly given by 


ftO= 204 y 


= 1 = 
an Fpl (0;¢)°+0 "8:9; (86950) + sy caixa *;0;(O;:C0;x)]- (3-4) 
Performing field redefinition Ç = ¢ + f(¢) generates third-order terms from the second- 
order action (2.9) that cancel relevant surface terms containing ¢. As a result, we find a 
useful relation [44—46] 


SOK] + SOK] = SPE] + Spurl]. (3.5) 


In terms of ¢, the total action is simply given by the second-order action (2.9) and bulk 
interaction (3.2). 


—12- 


Three-point functions of ¢ can be written schematically as 


(Ck: (7)Gkeo (7) Cues (T)) = (Cr: (7) Giro (7) Oks (7)) + redefinition terms, (3.6) 


where the first term is generated by bulk interaction (3.2) and the second term is boundary 
contribution at 7 from field redefinition. Higher-order correction to the expectation value of 
an operator O(7) is calculated by the in-in perturbation theory 


(O(r)) = (|r exp ( f i ar Hil’) ) O(r) Ir exp (-: / F ar‘ Hy!) ) , (3.7) 


where T and T denote time and antitime ordering. In this case, the interaction Hamiltonian is 
simply Hin, = — f AParlpuk, where Lpui is the integrand of Spuk. For three-point functions, 
the operator is Cx, (T)Ckə(T)Ck(T). First-order expansion of (3.7) reads 


(O(r)) =2 l i in O (3.8) 
E M A A is defined as 
(Cra (7) Chea (7) Shes (7)) = (207)°5 (k1 + k2 + k3) (Cr: (T)Cka(T)Cka(T)}- (3.9) 
A bispectram has bcil type E item be cpresed a 
(Ck: (T0)Cka (T0)Sks (T0)} = iN Ges (70)? Geo (70)? + [Cka (T0)|" [Ck (o)l? 
+14: (70)I?1Ces(T0)17| ; (3.10) 
where fan ie a Goeincient Walled nonlinear parameter. 


3.1 Slow-roll approximation 


In the SR approximation, the first line in cubic self-interaction (3.2) has O(e?) coefficients, 
while the second line has O(e?) coefficients. At the end of inflation, when the relevant modes 
are in the superhorizon limit, the leading term in the field redefinition (3.4) is the first term, 
because the other terms decay outside the horizon. Therefore, the third-order action becomes 


À 1 ; . 

Sunk [6] = Mp / dt dx a” EZ + EAL)? — zica + O), (3-11) 
with C = È + e2¢?/4. With the in-in perturbation theory (3.8), we obtain bispectrum of ¢ as 
Ht 

= 3.12 
(Ck: (70) Oko (70) Ckes (T0))) O(e2) 32M4 et (I kaka)? (3.12) 
kek? 
xX |= 5 k? + Ejs `> kik? + Se 
i iAj aes 


To arrive at bispectrum of ¢, we have to include contribution from the field redefinition 


(Ck: (70) C2 (70) Ces (T0))) = (Ck (70) Ce» (70) Ses (70)})) 


Helm) x2 [lce (70)I?1Gea(70)1? + [Cka (T0)|*|Cks (T0)1? + |r: (T0)l° Iks (T0)1?] - (3.13) 
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However, assuming €3 as a constant, €2(79) diverges as To — 0. The explicit time dependence 


is given by 
Ti —€3x T 
€2(T0) = E2x (2) X Ey (1 — €3; log 2) . (3.14) 


Tx Tx 


At leading order of €2(79), the bispectrum of ¢ is 


H* 


= wl 
(kx (70) Chea (T0) ies (70))) O(e2) 32M4 e (hi hoks)3 (3.15) 
+ . k2 k2 
3 2, i>j ifj 
x | (€2 — €1) x > k? + €x 2 kik; + 8e1x = 


In the squeezed limit, where one of the wavenumbers, kı, is much smaller than the 
others, kı — 0, the bispectrum becomes 


Ht 


{Ges (70)Gien(70)6—tea 70)}) = (2er + dea) sora Aka? 


(3.16) 


In terms of the mode function and spectral tilt, the squeezed bispectrum can be written as 


(Gres (70) Ges (70)6—Hea (T0))) = — (ns (k2, 70) — 1)|Gix (70)? Ck (70), (3.17) 


which is called Maldacena’s theorem. Here, we have shown Maldacena’s theorem by explicit 
derivation. It can also be proven by a model-independent method by considering the evolution 
of the small-scale perturbation in a rescaled dS background due to the long-wavelength 
perturbation [47]. 

Another important point discussed in [44] is the bispectrum induced by the last term 
in (3.2) 


1 a 
Sbuik [G] D Mà fa dex a eeecc. (3.18) 


Although it is SR-suppressed compared to (3.11), such term induces bispectrum that evolves 
outside the horizon. 


TO H* 3 3 3 
(Cki (70) Cko (70) Oks (70))) ole) = (core log 2) SDN Erkka) (k? + k3 + k3), (3.19) 


where x denotes horizon crossing of the total wavenumber kı + ko + k3. It shows that bis- 
pectrum of Ç evolves outside the horizon due to log 7 dependence. However, for bispectrum 
of ¢, there is an additional contribution from the field redefinition (3.13). After summing 
with the second term in (3.14), the log To dependence cancels out. Therefore, bispectrum of 
¢ does not evolve outside the horizon. 


3.2 Ultraslow-roll inflation 


During USR period, €; becomes extremely small because it scales as a~®. Also, ¿2 = 0 because 
€2 is a constant. Therefore, the contribution of the bulk interaction to the bispectrum is 
negligible and the leading contribution comes from the field redefinition. The relevant terms 
in the field redefinition are [48] 


: Ç, (3.20) 


_ £242 
f(Q)= 76 tF 


=J = 


which contribute 


(Ges (70)Gea(70)Gea(70))) = 2 (2 + 3) [Ii (70)|?1Gua(T0)I? + lea (70) 1a (70)1? 


+|Ck,(T0)|7|Ceg (70) || - (3.21) 


Factor 3 comes from ¢,(19) = 3H¢,(7), which can be derived from the USR mode function 
(2.52). Hence, we can read the local non-Gaussianity as 


6 5 
iil =e INE oe (3.22) 


In squeezed limit, the bispectrum becomes 
(Ges (70) tea (70) 6—tea (70))) = 6] Ce, (70)? lka (70), (3.23) 


which clearly violates Maldacena’s theorem (3.17). Indeed, for a shift-symmetric action of the 
inflaton, Maldacena’s theorem is generalized in [49, 50]. The squeezed bispectrum (3.23) can 
be reproduced from a generalized Maldacena’s theorem after substituting the USR condition. 


3.3 Ultraslow-roll to slow-roll transition 


In case there is a transition of the second SR parameter, €&2 becomes large during the transi- 
tion. Therefore, the leading term in the third-order action is [35, 48, 51-53] 


Sbutk[¢] = ma f dr dĉz a aE eaeh E POR |. (3.24) 


Because we are interested in transition from USR to SR period, relation between Ç and € at 
the end of inflation is € = Å + e2¢7/4. However, because |€o(79)| < 1, contribution from field 
redefinition to the bispectrum is much smaller than the bulk interaction (3.24). Thus, the 
bispectrum is given by 


{Cia (70) Con (70) Gheg(70))) = —2M3 ia dr ex(r)€h(r)a2(r) (3.25) 
x Im [Cr (70) Ceo (70) Shes (70) Ch, (7) CE (7), (T) + perm] . 


In the next part, we will evaluate the time integral for sharp and smooth transition of €2(T). 
Similarly as in Sec. 2.4, we follow the parametrization of €9(7) in [35]. 


3.3.1 Sharp transition 


Consider evolution of the second SR parameter as step function, given by (2.70), its time 
derivative is Dirac-delta function 


€,(7) = Acô (T — Ts). (3.26) 
Substituting it to the time integral in (3.25) leads to 


(Ger (0) Cho (70) Sieg (T0))) = — 2MP}€1 (Ts )a? (Ts) Aco (3.27) 
x Im [es (70) Cko (70) Ck (70) Ce, (Ts )Cka (Ts) Ch, (Ts) T perm] : 
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Then, we can substitute ¢,(7) and Ç; (T) given in Sec. 2.4.2 to obtain the bispectrum. Al- 
though the general form of bispectrum is complicated, in the limit k1, k2, k3 < ks the bis- 
pectrum becomes 


3 
(Chey (70) Gea (70) Cres (70))) = 2MPie1 (Ts)? (Ts) Ae (amas) 


< (eg gata) (3.28) 
k2 (KBK KBR RRJ’ i 


which is a local type bispectrum. Thus, we can read the nonlinear parameter fni = 5/8. 


3.3.2 Smooth transition 


Consider evolution of the second SR parameter given by (2.65), its time derivative is 


itis —2s?(s — 3)(s + 3) Ts\8 
(7) z o (3.29) 


T 


Because the second SR parameter is constant for T < Ts, the domain of time integral in (3.25) 
is simplified to 


litoe Gh (70))) = —2M3 [ dr e1(r)€5(r)a2(r) (3.30) 
x Im [Cr (70) Cho (T0) Cs (70) Ce, (7) Cia (7) ie, (T) + perm] . 


Unfortunately, the general form of bispectrum is complicated. However, it becomes simple in 
the squeezed limit kı < ko = kg for perturbation modes that cross the horizon much before 
the transition ki,k2,k3 < ks. Substituting (2.64), (3.29), and (2.68) to the time integral 
leads to 


(Chey (70) Chea (70) Cea (70))) = |r: (70)? lk (70) 7 [a Ts (r) +5200) 
w Peal TATE (3.31) 


Thus, we can read off the nonlinear parameter 


fnt = —2nv aay (3.32) 


which is suppressed by SR parameter ny. This means that the bispectrum generated by the 
transition (2.65) is negligible. In [54], it is argued that suppression in the bispectrum is a 
result of (2.65) satisfying Wands duality condition [55], which makes the integrand of (3.24) 
a total time derivative. The squeezed bispectrum in both smooth and sharp transition cases 
can be reproduced from the generalized Maldacena’s theorem by considering only background 
evolution [56] following [47]. 
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4 Generation of large fluctuation 


In Sec. 2.2, we have shown that violation of the SR approximation, especially the USR con- 
dition, leads to the growth of curvature perturbation outside the horizon [23, 24]. Although 
the USR period cannot realize the required e-folds number of inflation, it can be useful for 
other purposes. If the inflaton undergoes USR period only for a short duration, the cur- 
vature perturbation will be enhanced only within that duration. At the end of inflation, it 
corresponds to amplification of the curvature perturbation at a range of wavenumbers that 
cross the horizon during the USR period. 

Assuming sharp transitions, evolution of the second SR parameter corresponding to 
the potential in Fig. 1 is shown by the brown curve in Fig. 2. The inflaton undergoes USR 
period in the flat region of potential from time ts to te. In order to explain the observed CMB 
fluctuations, at early time, t < ts, the inflation was in SR period with slightly tilted potential. 
At late time, t = te, there is no observational constraint on the potential. A potential is viable 
as long as the curvature perturbation does not grow during the final period. 

In this section, we study the power spectrum and bispectrum of curvature perturbation 
in such a setup. 


4.1 Power spectrum 


At early time, t S ts, the inflaton was in SR period with Bunch-Davies initial vacuum, so 
the curvature perturbation is given by 


= iH —ikT . 
G(r) = (aca ma) $ (1+ ikr), (4.1) 


where €,, is the first SR parameter during this SR period. 

At ts St S te, the inflaton is in the intermediate USR. period. We define T, and Te as 
the conformal time corresponding to ts and te, respectively. The first SR parameter can be 
written as €1(T) = €15(T/7;)° based on proportionality in (2.47). Therefore, the curvature 
perturbation becomes 


E iH Ts 3 —ikT - ikT - 
G(r) = (zz =) © | Ave" (1+ ikr) — Bre” (1 — ikr)|, (4.2) 


where coefficients A; and Bp are determined by matching to the SR solution (4.1) at the 
boundary. Solutions of the coefficients by requiring continuity of ¢,(7) and (T) at transition 
T = Ts are 


3(1 + k?7?) 


3(1 + ikTs)? ikr, 


At late time, t > te, the inflaton goes back to SR dynamics. The curvature perturbation 
can be written as 


= iH Ts i —ikT . ikr z 
G(r) = (zz =) (2) [Cre ™ (1+ ikr) — Dyel*™(1—ikr)], (4.5) 


a 


where coefficients C and Dk are determined by matching to the USR solution (4.2) at the 
boundary. Solutions of the coefficients by requiring continuity of ¢,(7) and (T) at transition 
T = Te are 


— =l LaNa \2 2ik(Te—Ts ) 
Ch = IRI {9(krs i) (kTe + i)“e 
— [k?72(Qkrs + 3i) + 3i] [k?72(QkTe — 34) — 3i]}, (4.6) 
De = p gg (eB + h?r = Zikr) (er =i)? 
4kôrèrs 5 ° i 


+ie~2*re [3i + k?r2(2kTs + 3i)] (kre — it (4.7) 


As before, ks and ke are wavenumbers which cross the horizon at 7, and Te, respectively. At 
the end of inflation, To (— 0), the tree-level power spectrum is 


Aja) (ks 70) = = Gu (70)? = as i “ie — Dal? (4.8) 
a(0) rr Oe ggZ RNA 8r? M? e1 k k Bi“ f 
p s y S 


On large scale, the power spectrum approaches an almost scale-invariant limit 


A? ap (k) = A? (k < ks, To) = o 1+0 a (4.9) 
Aon an l 87? M2 Ers k2} |? 
* 


which is consistent with CMB observation. Assuming ke > ks, the next order approximation 


is 
H? 8 [ke\’ (k\? 24 (hke\? (kk \* kê 
Aso) (#70) = 1 ‘ + O . 
s(0) (E> 70) (r). | 5 (#2) (<) + 35 (F: ks k6 
(4.10) 
The minus sign in the coefficient of (k/k)? is important because it can explain a dip in the 
power spectrum before amplification start. As k approaches ks, the k?-term grows until it 


reaches unity. At this point, the k*-term is not large enough yet, so the power spectrum 
almost vanishes. This yields a dip in the power spectrum, which happens at the wavenumber 


1/2 
SeA 


It has been argued that this dip may lead to an interesting observational consequence [57]. 

At kaip < k < ks, the k*-term grows while the k®°-term is still small, which explains the 
apparent k4-growth in the power spectrum [38]. On a small scale with a larger wavenumber, 
ks Sk Ske, the power spectrum is oscillating around 


ke\ f 


whose high-density peak may collapse into PBHs. Numerical coefficient 4 appears for the 
same reason as (2.78). It is amplified by a factor of (ke/ks)f compared to the CMB-scale 


— 18 — 


ASPEN) i 


s0) (K) 


AŠ(sR) (px) 


Dx kaip ks ke 


Figure 3. Power spectrum of the curvature perturbation. At CMB pivot scale, p, = 0.05 Mpc™}, 


the power spectrum is almost scale invariant with amplitude A? (sr) (p+) = 2.1 x 1079, based on 
observational result [21]. At kaip < k < ks, the power spectrum grows with function k* until it 
reaches the peak. At ks < k < ke, the power spectrum is oscillating around 4 x AŽ PBN): Similar 
behavior at k > ke because SR period is imposed after the USR. period. 


power spectrum. To generate a significant abundance of PBHs, typically peak with © (0.01) 
on small-scale power spectrum is required. It means that the power spectrum is amplified 
O(107) compared to CMB-scale. In terms of the number of e-folds during the USR. period 
Nusr = log|a(te)/a(ts)], we need to have 


ke 1 
Nea log = = log 7“ ~ = log 10” = 2.68, (4.13) 
e S 


to achieve a significant abundance of PBHs. 


In this model, we can see that the power spectrum remains amplified even for k > ke, 
which oscillates around A? PBH)’ This happens because the inflaton goes directly to the SR 
period after the USR period. In case the inflaton undergoes a constant-roll period with e2 > 0 
after the USR period, the power spectrum at the end of inflation will decrease for k > ke. 


In this review, we choose an USR period for the intermediate phase as an example. In 
principle, the intermediate phase can be generalized to a constant-roll period as long as the 
curvature perturbation grows outside the horizon. The requirement was discussed below Eq. 
(2.63). However, for a general constant-roll period, the mode function becomes the Hankel 
function, which is given by (2.24) with (2.57). Technically, it becomes more difficult to 
match the coefficient of constant-roll period with the SR period because it involves derivative 
of Hankel function. In [40], a method to approximate the mode function is proposed by 
utilizing the fact that the duration of the SR violating regime is extremely short compared 
to the total duration of inflation. An inflation model with two phases of SR violating period 
is also introduced in [40] that yields k® growth in the power spectrum, which is much steeper 
than k* growth. 
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Figure 4. Plot of Co(k) and ns(k, 79) — 1. We choose ke/ks = 10 only for illustrative purposes. 


4.2 Bispectrum 


Consider the evolution of the second SR parameter shown by the brown curve in Fig. 2. Its 
time derivative can be modelled as Dirac-delta function 


€5(T) = Aeg[—6(7 — Ts) + 6(T — Te)]. (4.14) 


Substituting it to the time integral in (3.25) leads to 


(Cher (70) Geo (T0) Cie (70))) = 2Myje1 (Ts)? (Ts) Ae21m [Ck (7) Ca (7) Cea (70) Cie (Ts) Gio (Ts) Cts (Ts) 
—2Mpie1(Te)a7 (Te) AeoTm (Cx, (70) Ces (70) Chg (70) GE, (Te) City (Te) Ces (Te)] + perm. (4.15) 


In squeezed limit, the bispectrum becomes 


2 Ti 
tiie ata =< {4aeafe (eri An, 0) G rl) 
[Cka (T0)| 


Cin (T0) 


—4Ae3M 4e (Ts)a?(Ts)Im E c» | lr: (70) 171 Ceo (70) *- (4.16) 
[Cka (70) | 


We define the term inside parenthesis as Co(k), so the bispectrum can be written as 


(Cie, (70) Chea (70) Cea (70))) = —Co(k2)|¢k: (70)|7|Geo (70). (4.17) 


We can show that Co(k) = ns(k,7) — 1 analytically, which confirms Maldacena’s theorem 
explicitly [10]. 

At the end of inflation, we have shown that only the bulk interaction (3.24) contributes to 
the bispectrum. During the USR period, the bispectrum also satisfies Maldacena’s theorem, 
where both the bulk interaction (3.24) and field redefinition (3.20) are important. Therefore, 
in general time, the squeezed bispectrum is 


Keka (T)Cka(T)C-ka(T)} = — (ns (k2, T) — 1) (T)? lk (7) (4.18) 


We can read that spectral tilt of the short-scale perturbation, depends on time even at 
superhorizon scale, specifically around the transition time. Plot of the spectral tilt at the 
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Figure 5. Plot of n, — 1 at the end of USR. period and at the end of inflation. For a range 
of wavenumber between p, and ks, the spectral tilt evolves, altough such perturbations already far 
outside the horizon at the end of USR period. 


end of USR period and at the end of inflation are shown in Fig. 5. Thus, in contrast 
to the case of the SR inflation in Sec. 3.1 where time dependence log To is cancelled, the 
squeezed bispectrum evolves outside the horizon in this model. In case all wavenumbers have 
magnitude around the peak scale, the bispectrum is expected to deviate significantly from 
Maldacena’s limit. 

In the presence of peaks and troughs in the power spectrum, Maldacena’s theorem is 
generalized in [58]. In addition, the consistency relation (4.18) fixes the sign of fyn at the 
peak, which has an implication for the abundances of PBH [59]. 


4.3 Beyond tree-level contribution 


So far, we have discussed two-point and three-point functions of the curvature perturba- 
tion at tree-level approximation. We have shown that the sharp transition of the second 
SR parameter generates a O(1) contribution to the nonlinear parameter fng. It is a large 
contribution to the bispectrum, which implies that a higher-order correction than first-order 
perturbation theory can also be large. 

The second-order expansion of in-in perturbation theory reads 


(OTa = (O(7))fo.2) + OTa + (O(7))(0.2): (4.19) 
) 


(0,2 
(O(7)) ayy = [. dry f. dre (Hine (71) O(7) Hine (72) ) ; 
(O(T)) (0,2) =- I dr, i dry (O(7) Hine (71) Hine (72) ) l 


In the presence of cubic self-interaction in the Hamiltonian, the second-order perturbation 
contributes to one-loop correction to the power spectrum and tree-level contribution to the 
trispectrum. For an inflation model with sharp transition of the second SR parameter, a 
one-loop correction to the large-scale power spectrum was first derived in [9]. In this case, 
the operator is Cp(T)Ç-p(T), where p is the CMB scale wavevector, evaluated at T = To (— 0). 


= 9] = 


At the end of inflation, there is no contribution from the field redefinition because it 
is evaluated at the SR period. Only the bulk interaction (3.24) contributes to the one-loop 
correction. Substituting the corresponding Hamiltonian of (3.24) to (4.19) yields 


(Gp(70)6-»(70) aay = Må [. "nejat / "ertha 


«JI Re Jo (Fey + kz + ks)ð(ka + ks + ke) 
x (Chey (71) Chea (71) Ges (71) Sp (T0)6—p(T0) Shey (72) Ces (72) Ghee (T2))» (4-20) 


(Gp(70)6-»(70) 10. =- FMA S Tand aiai / daw Ghana) 


JI rs | oe + a + OCR + ks + ko) 
x (Cp(70)6—p(T0) Chey (71) Chea (71) Ges (71) Chey (T2)Cks (72) Chg (72))- (4-21) 


Performing Wick contraction and summing both terms lead to 


dèk 1 rls * * 
COl) = JECA f Se ONG GEGER 
+86) Colaba CnC + 86,6360 UEG + OGUE E)r —RelGp(r0)Go(7) (4.22) 
x (AGC Cala CG + GUGA + BOGE UG + OGUR E)n] Y 


Here, we only show the leading contribution to the one-loop correction that comes from 
integral of Dirac-delta function at T = Te. The contributions of T = 7, are subleading, which 
are explicitly evaluated in [10]. 

For p < k, the one-loop correction can be simplified to 


(olro) Gpl) Ya) = GMA (Aer) llo) lo 6 [eal P l] 


"(4.23) 
We find that the one-loop correction is exactly proportional to the quantum commutator 
(2.46), which implies the importance of decaying mode even for superhorizon mode function. 
If we only consider the constant mode of curvature perturbation at superhorizon scale, we 


will obtain an incorrect result [o(r), Ga] = 0. From commutator (2.46), we find 


1 


m(¢x(7)¢x' (T)) = IMZ) (4.24) 


which does not depend on wavenumber k. As we have already seen, in the standard SR 
inflation this quantity decreases exponentially as a function of conformal time. However, 
during USR regime, since € decreases much more rapicly than a~?, the right-hand-side of 
(2.46) and (4.24) increases exponentially, which makes the one-loop correction to the power 
spectrum large even on the large scale probed by CMB. 


= 9) = 


Then, substituting (4.24) into (4.23), we obtain a one-loop correction to the power 
spectrum 


1 

A3) (P, To) a gAe2)As(o)p, n) f| Ao (Es Te) (4.25) 
IR 

Because we are interested in the effect of large fluctuation on small scales related to PBH 

formation to the fluctuations on large scales observed by CMB, a natural choice of cutoffs are 

kir = ks and kyy = ke. For ke/ks < 1, evaluating the wavenumber integral with small-scale 

mode function (4.2) yields 


1 k 
A$) (P, To) = 1 (Ae2)’A2o) (p, 7) AS (Pem) (u + log r) : (4.26) 
Requiring the one-loop correction to be much smaller than tree-level contribution yields an 
upper-bound on small-scale power spectrum 


AŠ(PBH) < ~ 0.03, (4.27) 


1 
(Ae2)? 
where the numerical value is obtained by substituting Ae: = 6 for transition from USR to 
SR period. Such a bound suggests that sharp transition of the second SR parameter cannot 
generate large fluctuation on small scales, otherwise perturbation theory breaks down. This 
result should be understood under caveat that there are other contributions from quartic 
Hamiltonian, both induced by the third-order action (3.24) and fourth-order action. 


5 Discussion and summary 


5.1 Discussion 


Following [41], PBH formation models in single-field inflation can be classified into two cate- 
gories. The first category is models with features in the inflationary potential or nonminimally 
coupled inflaton with potential defined in the Einstein frame. Models with an extremely flat 
feature or an inflection point in the potential fall in this category [33, 60-83]. These mod- 
els can be derived from high-energy theories such as supergravity [60-66], axion monodromy 
[67, 68], scalaron in R?-gravity [69], a-attractor [70-73], and string theory [74-76]. Also, there 
are models from Higgs inflation which do not require theories beyond the standard model 
[77-81]. Other examples of features are a tiny bump or dip [84-89], an upward or downward 
step [32, 90-95], polynomial shape [96-99], and Coleman-Weinberg potential [24, 100-102]. 
In these examples, modification of potential makes a sharp transition of the second SR. pa- 
rameter on the inflaton dynamics. 

The second category is models with modified gravity or beyond nonminimally coupled 
inflaton. Examples for this category are models based on k [103] or G [104] inflation [105-110], 
the effective field theory of inflation [111-113], f(R) gravity [114], a nonminimal derivative 
coupling [115-117], Gauss-Bonnet inflation [118, 119], and bumpy axion inflation [120]. In 
these examples, amplification of small-scale perturbation can be realized by a sharp transition 
of the second SR parameter and/or other parameters. In [111, 112], amplification of small- 
scale perturbation is caused by a sharp transition of the sound speed, cs, a quantity that 
parametrizes deviation from a canonical kinetic term. In this case, coupling e% in cubic self- 
interaction (3.24) is modified to (€2/c2)! [121, 122]. For this category, rapid slowdown of the 
inflaton’s velocity might not be the only reason of amplification of the curvature perturbation. 
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There have been some works following [9, 10], most of them derive and qualitatively 
confirm the one-loop correction by other methods, although quantitatively has not been 
settled yet in the community. Despite of that, our claim [9] was criticized by the author of 
[123] who attempted to evaluate the same correction with a different method. In case of sharp 
transition of the second SR parameter, we have shown in [10] that after more careful and 
proper calculation, one can reproduce our result even if the method proposed in [123] is used. 
In case the second SR parameter evolves continuously with function (3.29), the bispectrum 
is suppressed by a SR parameter ny [124]. In [54], we have argued that suppression on the 
bispectrum and one-loop correction occur only for a specific smooth transition that satisfies 
Wands duality condition [55]. 

Indeed, there are issues related to regularization and renormalization of the divergence. 
In [9], we show that constraint (4.27) corresponds to the requirement to perform the renor- 
malization order by order in perturbation theory. Another method for performing regulariza- 
tion, renormalization, and resummation of loop corrections is proposed in [125-127], where 
the sound speed of the curvature perturbation is introduced within the effective field theory 
of inflation framework. Numerical integration of the one-loop correction has been performed 
in [128, 129], confirming our analytical result in the limit of the sharp transition. Also, Ref. 
[128] pointed out that a dip in the power spectrum (4.11) is a tree-level artifact, which disap- 
pears after adding a one-loop correction. One-loop correction to the power spectrum around 
the peak scale has been investigated in [130]. Specifically, for an inflationary potential with 
oscillatory features, it has been calculated in [131, 132]. 

Moreover, the 6N formalism and stochastic inflation in the presence of a sharp transition 
of the second SR parameter have been investigated in [133, 134]. The one-loop correction 
has also been derived from the ôN formalism in [135, 136]. The extension from USR to 
constrant-roll period has been derived in [137], which yields a constraint involving a constant 
second SR parameter and small-scale power spectrum. In the limit |€2| > 00, bispectrum and 
one-loop correction has been investigated in [138, 139], based on a method that approximates 
the curvature perturbation in the presence of a very short period where the SR approximation 
is violated [40]. In addition to cubic self-interaction, quartic self-interaction [140, 141] in the 
Hamiltonian can also contribute to the one-loop correction, in first-order perturbation theory. 
Its contribution has been investigated, either by direct expansion of the quartic derivative of 
the potential [142] or by an effective field theory approach [143, 144]. 

However, there are two works that completely disagree with our conclusions [145, 146]. 
In their approach, instead of performing field redefinition, they calculate a one-loop correction 
directly from third-order action with total time derivative terms and terms proportional to 
the equation of motion, as shown in (3.1). Currently, they have been criticized by [147]. 
Other perspectives on the effect of the total time derivative on the correlation function have 
been presented in [148, 149]. We hope to solve this discrepancy by pointing out what is 
missing in [145, 146]. Beyond the power spectrum, one-loop correction to the bispectrum has 
been explored in [150]. 

While we have been performing computation in the comoving gauge, some authors do 
so adopting the flat-slicing gauge. In this gauge, perturbation is represented by the inflaton’s 
fluctuation 6¢@, with vanishing curvature perturbation. Prior to our work, the quantum 
correction to 6d was investigated in [151, 152] and followed up in [153], which is consistent 
with our claim. The effect of a smooth transition of the second SR parameter in the flat-slicing 
gauge has been explored in [154, 155]. However, there is a result [156] that is inconsistent 
with our claim due to inappropriate order counts. 
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5.2 Summary 


In this chapter, we have explained the mechanism to generate large primordial fluctuation 
in single-field inflation. In order to amplify the curvature perturbation, violation of the SR 
approximation is needed. We have discussed the simplest example of an inflation model that 
violates the SR approximation, called USR inflation. We have shown that USR inflation 
leads to growth of the curvature perturbation outside the horizon. The consequence is that 
USR inflation cannot explain the required e-folds number to solve the horizon problem. 

Nevertheless, one can utilize the properties of the USR period for another purpose. If 
the USR period is sandwiched by two SR periods, the curvature perturbation grows tempo- 
rary during the USR period. At the end of inflation, the power spectrum of the curvature 
perturbation is amplified on small scales corresponding to the time scale of the USR period, 
while keeping the small amplitude on large scale consistent with CMB observations. 

Moreover, we have explained the nonlinear evolution of the curvature perturbation 
induced by cubic self-interaction, which yields three-point functions or bispectrum. The 
cubic self-interaction consists of bulk interaction and field redefinition. In USR inflation, 
only the field redefinition contributes to the bispectrum. When a transition from USR to SR 
period exists, the bulk interaction that is proportional to the time derivative of the second 
SR parameter becomes important. If the transition is sharp, the nonlinear parameter of 
the bispectrum has a O(1) value. This means that the bispectrum is large enough, which 
indicates that the nonlinear effect cannot be neglected. At the same time, such a nonlinear 
effect gives backreaction to the power spectrum as a one-loop correction that leads to a 
constraint on the small-scale power spectrum associated with transition of the second SR 
parameter. 
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